First catalog of strong lens candidates in the COSMOS field 



Cecile Faure 1 ' 2 , Jean-Paul Kneib 3 ' 4 , Giovanni Covone 3 ' 5 , Lidia Tasca 3 , Alexie Leauthaud 3 , 
Peter Capak 4 , Knud Jahnke 11 , Vernesa Smolcic 11 , Sylvain de la Torre 3 , Richard Ellis 4 , 
Alexis Finoguenov 10 , Anton Koekemoer 7 , Oliver Le Fevre 3 , Richard Massey 4 , Yannick 
Mellier 8 , Alexandre Refregier 9 , Jason Rhodes 4 , Nick Scoville 4 , Eva Schinnerer 11 , James 
Taylor 4 ' 12 , Ludovic Van Waerbeke 6 , Jakob Walcher 3 



ABSTRACT 

We present the first catalog of 67 strong galaxy-galaxy lens candidates discovered in the 1.64 
square degree Hubble Space Telescope COSMOS survey. Twenty of these systems display multiple 
images or strongly curved large arcs. Our initial search is performed by visual inspection of the 
data and is restricted, for practical considerations, to massive early-type lens galaxies with arcs 
found at radii smaller than ~5". Simple mass models are constructed for the best lens candidates 
and our results are compared to the strong lensing catalogs of the SLACS survey and the CAS- 
TLES database. These new strong galaxy-galaxy lensing systems constitute a valuable sample 
to study the mass distribution of early-type galaxies and their associated dark matter halos. We 
further expect this sample to play an important role in the testing of software algorithms de- 
signed to automatically search for strong gravitational lenses. From our analysis a robust lower 
limit is derived for the expected occurrence of strong galaxy-galaxy systems in current and fu- 
ture space-based wide-field imaging surveys. We expect that such surveys should uncover a large 
number of strong lensing systems (more than 10 systems per square degree), which will allow for 
a detailed statistical analysis of galaxy properties and will likely lead to constraints on models of 
gravitational structure formation and cosmology. The sample of strong lenses is available at this 
address: |http: / /cosmosstronglensing.uni-hd.de/| 

Subject headings: catalogs — gravitational lensing — galaxies: statistics 
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1. Introduction 

Gravitational lensing is one of the most promis- 
ing tools in modern astrophysics for probing 
the matter content of the universe directly, and 
thereby constraining cosmological models. Ap- 
plications of strong gravitational lensing include 
measurements of the Hubble constant using lensed 
quasars with time-delays (Refsdal 1964; e.g. of 
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time delay measurements: Barkana 1997; Co- 
hen et al. 2000; Burud et al. 2000; Burud et 
al. 2002 a,b; Fassnacht et al. 2002; Hjorth et 
al. 2002, Ofek & Maoz 2003; Jakobsson et al. 2005; 
Kochanek et al. 2006; Fohlmeister et al. 2006; 
Vuissoz et al. 2007) and the determination of 
galaxy and cluster mass distributions using mul- 
tiple images (e.g. Guzik & Seljak 2002, Mandel- 
baum et al. 2006, Read et al. 2007). Moreover, 
homogeneous samples of strong lenses allow for a 
statistical determination of the properties of lens- 
ing galaxies (e.g. Kochanek 1996; Rusin & Ma 
2001; Chae, Mao & Kang 2006, Chae 2007) and 
samples of strong lenses selected based on source 
properties are used to test cosmological models 
(e.g. Fukugita & Turner 1991; Wambsganss et 
al. 2005; Moeller, Kitzbichler & Natarajan 2006; 
Hilbert et al. 2007). 

The giant arcs that are observed around mas- 
sive clusters of galaxies are among the most spec- 
tacular lensing phenomena (e.g. Gladders et 
al. 2002 searched giant arcs in the Red-Sequence 
Cluster Survey (RCS); Scarpine et al. 2006, Hcn- 
nawi et al. 2006, Estrada et al. 2007 in the Sloan 
Digital Survey (SDSS)). These cases are relatively 
rare because each square degree of sky contains 
only about one cluster of sufficient mass to pro- 
duce a giant arc. Much more common are the 
strong lensing events found around massive early- 
type galaxies. The first studies of galaxy-galaxy 
strong lensing systems targeted distant, bright 
lenscd sources, including quasars (ESO Hamburg 
survey: Wisotzki et al. 1996; CfA-Arizona-ST- 
LEns-Survey: Peng et al. 1997, Muhoz et al. 1998) 
and radio sources (Jodrell Bank-VLA Astrometric 
survey: Patnaik et al. 1992, Browne et al. 1998; 
Cosmic Lens All Sky Survey: Jackson et al. 1995, 
Myers et al. 1995). The identification of multiple 
images in these systems is generally unambiguous 
and does not require very deep observations; but 
requires imaging over a large fraction of the sky. 

The advent of deep and wide optical imag- 
ing combined with large spectroscopic surveys has 
greatly broadened the opportunities and means by 
which strong galaxy-galaxy lensing events can be 
found. Indeed, since we expect that ~0.1% of 
massive galaxies (the luminous elliptical galaxies) 
strongly distort and magnify a galaxy source at 
z>l (Miralda-Escude & Lehar 1992), new wide 
field surveys are likely to discover many hun- 



dreds of fainter galaxy-galaxy strong lensing sys- 
tems. Spectroscopic identification of two galax- 
ies at different redshifts in the same Sloan Digi- 
tal Sky Survey fiber has proved to be a successful 
method in the Sloan Lens ACS survey (SLACS: 
Bolton et al. 2006). Deep, wide-field imaging sur- 
veys (especially those utilizing the exquisite im- 
age resolution available from space) provide an- 
other avenue for locating strong galaxy-galaxy 
lenses. This was first attempted in the Hubble 
Space Telescope Medium Deep Survey (Griffith et 
al. 1994, Ratnatunga et al. 1999), and more re- 
cently in the Great Observatories Origins Deep 
Survey field (GOODS; Dickinson et al. 2001, Fass- 
nacht et al. 2004) or in the AEGIS survey (Mous- 
takas et al. 2006). Automated software for identi- 
fying strong lenses in HST image archives are cur- 
rently under development (Marshall et al. 2005, 
Moustakas et al. 2006) and appear to be very 
promising. Searches for strong lens systems are 
also underway from ground-based wide-field sur- 
veys. In particular the SL2S initiative (Cabanac 
et al. 2007) was able to find a few lenses per square 
degree in the CFHT-LS wide survey. However, 
ground based imaging surveys are necessarily lim- 
ited to finding systems with large deflection an- 
gles because of the seeing limitation compared to 
space-based observations. 

The 1.64 square degree Hubble Space Telescope 
COSMOS survey (Scoville et al. 2007) provides an 
excellent opportunity to locate and study a large 
number of strong galaxy-galaxy lensing systems. 
It includes the largest contiguous high-resolution 
astronomical imaging survey ever performed from 
space. Although observations with the HST Ad- 
vanced Camera for Surveys (ACS) have only been 
obtained in a single band (F814W filter), multi- 
wavelength coverage of this equatorial field is pro- 
vided in the optical band by deep multicolor imag- 
ing from Suprimc on the Subaru telescope and 
Megacam at CFHT. 

In this paper, we present the first systematic at- 
tempt to locate strong galaxy-galaxy lensing sys- 
tems in the COSMOS field. The data-set used for 
our investigation is described in § [2] and the tech- 
niques used to uncover strong galaxy-galaxy lens- 
ing candidates are discussed in § [3] Our findings 
are presented in § [4] The lens potential modeling 
scheme is discussed in §[5] and we compare our lens 
sample to other strong lens catalogs in §[6] Finally, 
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the results are discussed in § [7J Throughout this 
paper, galaxy magnitudes are quoted in the AB 
system and cosmological distances are calculated 
assuming a ACDM model with parameters Ho = 
73 km s" 1 Mpc" 1 , ft m =0.3, and ft A =0.7. 



2. Data 



To conduct our strong galaxy-galaxy lens search 
we take advantage of three key data-sets of the 
COSMOS survey: the ACS/ HST high resolution 
imaging, the Suprime/SVt&arit imaging, and the 
Megacam/ CFHT multicolor imaging. 

2.1. HST data 

We have analyzed the full ACS field of view 
of the COSMOS field which corresponds to 1.64 
square degree observed during cycle 12 and 13. 
We use data release vl.2 of the COSMOS ACS 
data which can be accessed through the dedi- 
cated database at IRSAP] Individual exposures 
were combined using the MultiDrizzle software 
(Koekemoer et al. 2007) onto an output grid of 
pixel size 0.03". The limiting magnitude is about 
If814W = 26.5 mag. A detailed description of the 
data is available in Scoville et al. (2007) , and a de- 
tailed description of the ACS-selected galaxy cat- 
alog can be found in Leauthaud et al. (2007). 

2.2. Subaru and CFHT data 

Multicolor imaging of the COSMOS field 
was obtained from the Subaru telescope with 
the Suprime camera (Miyazaki et al. 2002) in 
Bj,Vj,g+, r+, i+, z+ bands, and from the CFHT 
with the Megacam camera (Boulade et al. 2003) 
in the u* and i* bands. In these data, the seeing 
is generally better than 1". A photometric cata- 
log was derived from a combination of these data, 
based on detections in the Subaru i+ band (Capak 
et al. 2007). 

2.3. Photometric redshifts 

The publicly available Le Phare photometric 
redshift estimation code (Ilbert et al. 2006) has 
been used to measure the redshifts of 278 526 
galaxies with Ip &uw < 25 mag. Details concern- 
ing the multi-wavelength photometry can be found 
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Fig. 1.— The galaxies of the COSMOS field: 
Vj-band absolute magnitude versus photometric 
redshift. The grey solid line contours show iso- 
density numbers of all the galaxies in the entire 
photometric catalog (levels: 1, 100, 300, 500, 1000, 
1200, 1600 galaxies). The black solid line contours 
draw the isodensity numbers of galaxies in the 'P'- 
catalog (levels: 1, 10, 30, 50, 100, 140 galaxies). 
The red crosses point the lensing galaxy candi- 
dates described in Sect. [4] 



1 http: / /irsa. ipac.caltech.edu/data/COSMOS/images/ 
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in Mobasher et al. (2007). The I F %uw < 25 mag 
limit ensures a good accuracy of both the photo- 
metric redshift estimation and the fit to the spec- 
tral energy distribution (SED). The ground-based 
photometric zero-points were calibrated (and the 
SED templates adjusted accordingly) using 1095 
spectroscopic redshifts from the zCOSMOS Sur- 
vey (Lilly et al. 2006). Using 8 bands, this method 
achieves a photometric redshift recovery accuracy 
of a Az /(l + z s ) = 0.031 with r] = 1.0% of catas- 
trophic errors, defined as Az/(1 + z s ) > 0.15 (fur- 
ther details can be found in Mobasher et al. 2007). 

3. Methodology 

A four-step procedure is employed in order to 
identify strong galaxy-galaxy lens systems. In or- 
der, the four steps are the following: 

1 . Select a list of potential lenses from the pho- 
tometric redshift catalog (hereafter 'P' cat- 
alojQ. 

2. Visually inspect the ACS images of all galax- 
ies in the 'P' catalog to produce a catalog of 
potential strong galaxy-galaxy lens systems 
(hereafter 'E' catalog |^] Only 302 galaxies 
are in this catalog, the 35 others are missing 
from our personal archives.). 

3. Investigate the 'E' catalog using multicolor 
images to check for color differences between 
the main galaxy and the potentially lcnscd 
object. 

4. Subtract a galaxy surface brightness model 
of the foreground galaxy to determine the 
morphology and lensing configuration of the 
potentially lensed background galaxy. 

In our methodology, systems that pass steps 
both 3 and 4 qualify as strong galaxy-galaxy lens 
candidates. We will now describe each step in 
more detail. 

3.1. The potential lens catalog 

The ACS-based photometric redshift catalog 
limited to F8UW = 25 mag contains 278 819 en- 
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2 This catalog can be found at 

heidelberg.de/mitarbeiter/cfaure/downloads.html/Pcatalcgxat] 
3 This catalog can be found at |http://www.a ri.uni- 



heidelberg.de/mitarbeiter/cfaure/downloads.html/Ecatalog.cat 



tries: it is unrealistic to inspect visually all galax- 
ies in this catalog. Furthermore, this would be 
inefficient, because we know that a galaxy needs 
to be massive and not too distant to efficiently dis- 
tort and strongly magnify a more distant galaxy. 
We therefore limit our search to intrinsically bright 
galaxies at moderate redshifts. Furthermore, for 
this first investigation of the COSMOS field, we 
only look at early-type galaxies (classified spec- 
troscopically by the photometric redshift code), 
which generally have a simple, smooth surface 
brightness distribution. Once again, we empha- 
size that the main reason for applying these selec- 
tion criteria to the photometric redshift catalog is 
to reduce the sample of galaxies that need to be 
visually inspected. It is highly possible that ad- 
ditional strong lens systems will be found around 
galaxies that do not belong to our search sample. 
The number of strong lensing systems that we find 
is thus a lower limit to the total number of strong 
lensing systems in the COSMOS field. 

To search for potential lenses, we select galaxies 
according to the following criteria: 

• Their photometric redshift must be in the 
range 0.2 < z p hot < 1-0. The lower redshift 
bound of z = 0.2 is chosen to ensure that 
galaxies will indeed produce lensing events 
and the upper bound is essentially motivated 
to ensure a good accuracy on the photomet- 
ric redshift given that no deep infrared data 
are yet available for the COSMOS field. The 
galaxy is retained if the error bars at 68% 
confidence level overlap, even partially, the 
redshift range [0.2; 1.0]. 

• Their luminosity must be My < —20 mag, 
on the assumption that the brightest galax- 
ies are also likely to be the most massive ones 
with the largest cross-section for lensing. 

• We limit our search to early-type galaxies 
that have been spectrally classified as such 
by the SED fitting of our photometric red- 
shift algorithm. The motivation for this cri- 
terion is to focus on massive systems with a 
simple light distribution that is easy to fit 
and to subtract. 

We obtain a final 'P' catalog of 9452 bright el- 
liptical galaxies (see Fig. [TJ. Further strong lens- 
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ing searches using automated software will cer- 
tainly relax these criteria. 

During the course of our visual inspection, we 
serendipitously found some additional strong lens 
systems in which the lensing galaxy did not com- 
ply with the above limits. Although these con- 
stitute an inhomogeneous sample, we add them to 
the list of strong lensing systems for completeness. 
Those systems are identified in Tables [T] and [2] 

3.2. Visual inspection of the ACS images 

We inspect visually a postage-stamp 10"xl0" 
ACS image surrounding all galaxies in the 'P' cat- 
alog, looking for lensed features (tangential arcs or 
multiple images). The search box size limits us to 
T arc < 5". Consequently, this puts an upper limit 
on the mass scale for the lens, which excludes clus- 
ters as a possible main deflector. This inspection 
was performed by 5 of the co-authors using the 
DS9 FITS image visualization tool. 

Our strategy reduced the 9452 galaxies in the 
! P' catalog to an 'E' catalog of 337 candidates, 
with intentionally loose selection criteria to avoid 
excluding any real lenses at this stage. 

3.3. Pseudo-color images 

A color image is always useful to verify poten- 
tial lens candidates. Indeed, gravitational lensing 
is achromatic, thus a gravitational arc or multi- 
ple images systems should have a consistent color 
across the image(s). Furthermore, the color of a 
lensed galaxy will only rarely have the same color 
as the main lensing galaxy. This color test hence 
provides a good way to discriminate false strong 
lensing cases. However, with our particular com- 
bination of ground- and space-based data (espe- 
cially when the seeing approaches 1"), the color 
test is only useful for arcs more distant than ~0.5" 
from the lensing galaxy, because of strong con- 
tamination from the wings of the lensing galaxy. 
For the 337 systems in the 'E' catalog, we pro- 
duced pseudo-color images made from the "best 
PSF" stacks of the B, r and i Subaru images, us- 
ing the RGB feature in the DS9 software, as well 
as ground-based images sharpened with the ACS 
images. The color images use the ACS F814W 
data as an illumination map and the Subaru Bj, 
r + , and z + images as a color map. To achieve 
this each Subaru image is divided by the average 



of the three Subaru images, then multiplied by the 
ACS F814W image. This preserves the flux ratio 
between images, while replacing the overall illumi- 
nation pattern with the F814W data. Each image 
is then smoothed by a 1 pixel Gaussian to reduce 
the noise and divided by A 2 to enhance the color 
difference between star forming and passive galax- 
ies. The processed Bj, r + , and z + images are 
then assigned to the blue, green, and red chan- 
nels respectively. The resulting images have the 
high spatial resolution of the ACS imaging, but 
color gradients at ground based resolution. Of the 
337 candidates, 67 have plausible colors for strong 
lensing systems. The color images are displayed 
in Figs. [3] to [13] For illustrative purposes, two 
objects that were rejected at this stage are shown 
in Fig. [2] The system displayed in the top panel of 
Fig. [2] shows two arc-like features in the ACS im- 
age: a giant arc located North- West of the main 
galaxy, and its candidate counterpart, located to 
the South- West. The pseudo-color image reveals 
that the "arcs" have very different colors, indicat- 
ing that they are not images of the same source: 
neither arc likely to be distorted images of back- 
ground galaxies. Pronounced, arc-like features are 
also observed in the ACS image shown in the bot- 
tom panel of Fig. [2] North-East of the main galaxy 
bulge. With the single-color ACS data-set only, it 
would not be clear whether these features were 
part of the galaxy's spiral arms, or distorted im- 
ages of a background object. The pseudo-color im- 
age shows blue knots of star formation all around 
the disc of this galaxy, suggesting that the arc-like 
features are most probably genuine regions of the 
main galaxy. 

3.4. The lensing galaxy fit and subtraction 

Fitting and extracting the foreground lens 
galaxy light from the images reveals the shape 
of the more distant, lensed object. For this work, 
a robust 2-dimensional fit to the galaxy surface 
brightness distribution is derived using GIM2D 
(Simard 1998, Marleau & Simard 1998). We adopt 
a Sersic bulge plus exponential disc parametri- 
sation to describe the two-dimensional surface 
brightness distribution of the lensing galaxy light 
profile. The Sersic law is parametrised by means of 
the total flux in the bulge, the Sersic index n, the 
bulge ellipticity (e = 1 — b/a), the position angle 
of the bulge and the effective radius of the bulge 
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R e f f . The exponential profile depends on the pho- 
tometric disk total flux, the disk scale-length, the 
disk position angle and the disk inclination. The 
GIM2D software finds the best fitting values for 
all of these parameters, and provides them as an 
output of the decomposition process. 

In addition to the structural parameters for 
each object, we calculate a model image of each 
galaxy, which is convolved with a point spread 
function before comparison with real data. An im- 
age of the residuals is obtained by the subtraction 
of the model from the science image. We present 
the result of the galaxy fits in Tables [T] [2] and in 
Figs.[3lto[l2| Inmost cases, the arc candidates ap- 
pear clearly after the galaxy model is subtracted 
(panel 3 in those figures). 

In 10 cases, the modeling and subtraction of 
the lens galaxy proves difficult. The problems 
encountered include irregular galaxy morphology 
and compact galaxy size. We display the ACS 
and color images of these systems in Fig. 13 Ad- 



ditionally, the system COSMOS 5737+3424 was 
discovered serendipitously in a part of the ground- 
based data-set not covered by the ACS images. We 
present just its pseudo-color image in Fig. [13] 

4. Results of the strong lens search 

The location of lens candidates in the absolute 
V-band magnitude versus redshift plane are indi- 
cated by crosses in Fig. [T] 

We divide the sample in two sub-samples: the 
multiple arcs and long curved image systems, 
which we refer to hereafter as the "best systems" , 
and less certain single arclet systems. Given their 
distinctive shape, the multiple images and long 
curved arcs have a greater probability of being 
genuine strong lenses relative to the single ar- 
clet systems. Throughout the rest this paper we 
therefore restrict our study to the sample of "best 
systems" . 

The stamp images of the best systems are 
shown in Figs. [3j [4] [5] and [13] and discussed in 
Sect . |4.1| Images of the single arclet systems are 
displayed in Figs. [§|7][8j|9][lO][TT][l2]and[l3] 
and discussed in Sect. 14.21 Their characteristics 
are given in Tables [T] and [2] In Fig. [14] we have 
reported the distribution of the arc radius of the 
COSMOS strong lens candidates. Most systems 
have an arc radius< 2". 



4.1. The best systems 

In twenty cases, we have detected systems with 
a single long curved image or with two or more 
arcs of similar color symmetrically located around 
the lensing galaxy or with ring-like features. Their 
coordinates and characteristics are listed in Table 
[l] We comment on the most interesting or/ and 
less trivial systems hereafter. 

A bright arc-like feature is clearly visible in the 
ACS image of the system COSMOS 0012+2015. 
The counter-arc becomes visible after subtraction 
of the lensing galaxy profile, South-West of the 
center of the lens galaxy. 

The system COSMOS 0018+3845 is most prob- 
ably a triply or possibly quadruply imaged back- 
ground galaxy. The North-West arc would be a 
single image of the galaxy while the South-East 
arc would be a composite of two or three images 
of the same galaxy. 

The system COSMOS 0038+4133 is a complete 
Einstein ring candidate, formed by the distorted 
image of a single background galaxy. 

Four systems show ring-like structures around 
a central elliptical object (COSMOS 0049+5128, 
COSMOS 0124+5121, COSMOS 5941+3628 and 
COSMOS 5947+4752). In those cases, it is dif- 
ficult to conclude without any spectroscopic data 
whether the ring is due to gravitational lensing or 
if it is a genuine part of the central galaxy (like 
a ring-galaxy or a shell). The main argument in 
favor of them being lens candidates is the strong 
color difference between the elliptical central ob- 
jects and the ring features, suggesting they are not 
physically related. 

There are probably four images of a back- 
ground source in the candidate system COS- 
MOS 0050+4901. Three images are combined 
in an arc South-West of the lensing galaxy, while 
a counterpart image is visible North-East of the 
galaxy. 

The system COSMOS 0211+1139 displays two 
giant arcs, each with different colors in the Subaru 
image. They are probably distorted images of two 
different background galaxies. 

The pseudo-color image of COSMOS0216+2955 
shows a strong color difference between the galaxy 
and the potential arcs. Nevertheless, the galaxy 
subtraction shows features that are not fitted by 
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Fig. 3. — The first strong lensing candidates in the COSMOS field: the best systems. From left to right: 
the first panel is the ACS I-band image of the lens. A letter (f) close to the name of an object indicates 
that it was found fortuitously. The second panel is the Gim2D model of the lensing galaxy. The third panel 
is a subtraction of the two first images. The fourth panel is the color image of the systems from B,R,I 
ground based images. The fifth panel is the color image sharpened with the ACS-F814u; images. The label 
on the left top corner of panel five indicates the image multiplicity. When two multiplicities are given and 
are separated by a slash, the second multiplicity is possible but difficult to confirm with the present dataset. 
North is to the top of the image and East is to the left. The photometric redshift of the lensing galaxy is 
displayed in the second panel (error bars are given in Table [TJ and the radius of the arc is displayed in the 
third panel. 
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Fig. 6. — The first strong lensing candidates in the COSMOS field: the single arclet systems. From left 
to right: the first panel is the ACS Fband image of the lens. A letter (f) close to the name of an object 
indicates that it was found fortuitously. The second panel is the Gim2D model of the lensing galaxy. The 
third panel is the difference between the two first images. The fourth panel is the color image of the systems 
from B,R,I ground based images. The fifth panel is the color image sharpened with the ACS-F814w images. 
North is to the top of the image and East is to the left. The photometric redshift of the lensing galaxy is 
displayed in the second panel (error bars are given in Table |2| and the radius of the arc is displayed in the 
third panel. 
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Fig. 7. — Continue Fig. [6] 
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Fig. 9. — Continue Fig. [6] 
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Fig. 12.— Continue Fig. [6] 
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Fig. 13. — Other strong lensing candidates. For these systems, no satisfying fit of the lensing galaxy was 
possible. Therefore we display only the ACS I-band image of the system (left panel when available), the 
pseudo-color image (middle panel) and the color image sharpened with the ACS image (right panel). The 
labels refer to the multiplicity of the images when different from single and when known. The letter (f ) close 
to the name of an object means that it was found fortuitously. North is to the top of the image and East is 
to the left. 
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Fig. 14. — Distribution of the arc radii in the sample of COSMOS strong lens candidates (solid line his- 
togram). The striped histogram indicates the distribution of the best systems. 




6"x6" 



Fig. 2. — Examples of strong lens candidates 
ruled out by the pseudo-color image test. North 
is to the top of the figure and East is to the left. 
The left stamp images are the ACS I-band images 
of the systems. The right stamp images are the 
color images. 



the Sersic profile. They seem to be linked to the 
galaxy core. Hence, there are two possible ex- 
planations for this system: either it is actually 
a strong lens, and the lensing galaxy is a spiral 
galaxy, or the two arcs observed are part of the 
arms of a spiral galaxy. 

The system COSMOS0254+1430 has two lens- 
ing galaxies at very short distance. We have not 
been able to make a clean galaxy light profile fit. 
The two arclets located North-East and South- 
West have are clearly visible in the color image. 

In the system COSMOS 5737+3424, we observe 
at least three giant arc like features around the 
central galaxies of an identified XMM detected 
cluster (Finoguenov et al. 2007). 

COSMOS 5921+0638 is a quadruply-imaged 
lensed quasar or very compact source. An Einstein 
ring is clearly visible in the ACS/F814W band im- 
age. 

4.2. The single arclet systems 

Forty seven of the discovered lens candidates 
have a single arclet identified in both the ACS and 
the color images. In these cases, no counter image 
for these arcs could be found in the residuals of 
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Table 1 

Summary of the lensing galaxy parameters for the best systems 



Name 


RA 






DEC 






z±Az 


Tare 




n 


1' 


PA 


* 


5 


i 814tu 




h:m:s 




d:':" 








II 




II 










deg 










0012+2015 


10:00: 


12 


.6 


+02:20: 


:15 


0. 


4l ±aoi 


0.90 





.53 


8 


.1 


19. 


98-1-0.01 
^ 8± 0.01 


-10.2 


0. 


,25 


21 


.6 


0013+2249 


10:00: 


:13 


.9 


+02:22. 


:49 


0. 


oo_i_0.06 
■ oo - n 0.02 


1.65 





.90 


3 


.9 


18, 


■89±8:§i 


38.2 


0, 


16 


22 


.4 


0018+3845^ 


10:00: 


:18 


.4 


+02:38: 


:45 


0. 


■7l±§ :i | 


0.40 





.29 


5 


.6 


23, 


•6o±8:8i 


-21.8 


0, 


23 


22 


.7 


0038+4133 


10:00 


:38 


.2 


+02:41: 


:33 


0. 




0.74 





.72 


5 


.2 


20, 


oq+0.05 

oy=c 0.05 


-2.9 


0, 


26 


20 


.5 


0047+5023 


10:00 


:47 


.6 


+01:50 


:23 


0. 


CK+0.07 


0.70 





.72 


1 


.9 


20, 


GC + 0.01 

D0=t 0.01 


33.4 


0. 


,05 


22 


.8 


0049+5128 


10:00: 


:49 


.2 


+01:51: 


:28 


0. 


oo_i_0.03 
OOIC 0.05 


2.22 


1 


.07 


1 


.2 


19, 


R1+0.01 
D1=c 0.01 


-24.6 


0. 


22 


23 


.3 


0050+4901^ 


10:00 


:50 


.6 


+02:49 


:01 


1. 


01 +°- 07 

Ui=c 0.05 


1.90 





.36 


2 


.9 
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79-L0.02 
' z=c 0.03 


26.6 


0, 


39 


22 


.7 
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0. 
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4 


.6 
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0. 
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.3 
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.5 
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0. 
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2 


.4 
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0. 


,25 


23 
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11 


.2 
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0. 
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.8 
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22 
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.5 
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:51 
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6 


.3 


21, 


Q4+0. 03 
y4± 0.03 


80.9 


0, 


33 


22 


,3 


0254+1430 


10:02 


:54 


.0 


+02:14 


:30 






1.54 














-56.5 


0, 


,40 


22 


.1 


5737+3424^ 


09:57: 


:37 


.0 


+02:34: 


:24 






12.0 
























5758+1525^ 
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.6 


+02:15 


:25 






0.35 




















21 


.7 


5857+5949 


09:58: 


:57 


.0 


+01:59 


:49 


0. 


oqiO.03 
,3y=t 0.05 


2.15 





.61 


1 


.4 


20, 


nciO.Ol 
uo=c 0.01 


-32.8 


0, 


36 


21 


.9 


5914+1219 


09:59 


:14 


.7 


+02:12 


:19 


1. 


05+ - 05 


1.86 





.27 


3 


.9 


23, 


9CI0.80 

zo=c 0.04 


14.4 


0, 


,08 


22 


.8 


5921+0638 


09:59 


:21 


.7 


+02:06 


:38 


0. 


4C_i_0.03 


0.80 





.49 


6 


.8 


20, 


041O.O2 
,34:=C: 0.02 


29.9 


0, 


,08 


20 


.6 


5941+3628 


09:59 


AI 


.3 


+02:36 


:28 


0. 


qn+o.06 

yu=c 0.02 


1.21 





.76 


1 


.0 


20, 


Q1 +0.02 
yl=c 0.02 


-5.4 


0. 


,08 


22 


.8 


5947+4752 


09:59 


:47 


.8 


+02:47 


:52 


0. 


9010.04 


2.55 





.51 


1 


.6 


19, 


oo-lO.04 
■° J± 0.05 


4.1 


0. 


,07 


22 


.8 



Note. — Column 1: Lens candidate name. A letter * indicates that the object was found fortuitously. 
Columns 2 and 3: Coordinates in J2000. Column 4: Photometric redshift of the lens and error bars (at 
68% confidence level). Column 5: Largest radius of the arc in arc-seconds. Column 6: Effective radius of 
the lensing galaxy from the surface brightness fitting, in arcsec. Column 7: Sersic profile index. Column 
8: Total magnitude of the lensing galaxy in the WFPC/F814w band. Column 9: Lensing galaxy position 
angle. Column 10: Galaxy ellipticity e. Column 11: Magnitude of the brightest image in the WFPC/F814w 
band (in magnitudes per square arc-second). 
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Fig. 15. — The spectroscopic redshift distribu- 
tion of the lensing galaxies in the SLACS sur- 
vey (hashed histogram) and in the CASTLES 
database (solid line histogram), compared to the 
photometric redshift distribution of the best sys- 
tems in the COSMOS survey (filled-in grey his- 
togram) . 



the galaxy-subtracted image. As observed on the 
Subaru data, the arclets are generally much bluer 
than the lensing galaxies, and are likely to be only 
weakly distorted by the galaxy. Their coordinates 
and characteristics are give in Table [2] We com- 
ment below on a few interesting systems. 

The system COSMOS 0009+2455 shows two 
elongated features symmetrically distributed around 
the elliptical galaxy. The two features have 
slightly different colors, therefore if the system 
is a genuine lens, the sources are probably two 
different objects. 

The arclet in the system COSMOS 0028+1919 
is located West of the lensing galaxy, and appears 
in blue in the color image. 

In the system COSMOS 5805+0413, a giant 
blue arclet appears between two galaxies. The arc 
is slightly curved in direction of the western galaxy 
but there are no counterpart images around any 
of the two galaxies. The parameters of the early 
spectral type galaxy at the center of the stamp 
image are presented in the Table [2j 
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4.3. Radio and X-ray counterparts 

Six of the lens candidates have radio counter- 
parts in the 1.4GHz VLA-COSMOS data (Schin- 
nerer et al. 2006). The radio data cover the full 
COSMOS field to a depth of ~10.5/iJy/6eam rms 
in the center of the field. All six sources lie above 
a 5ct level in total flux, ranging from 77 /iJy to 
1.24 mJy in total 1.4GHz flux density. The dis- 
tance between optical and radio position is in all 
cases less than 0.5". 

One of the lens candidates, COSMOS 0013+2249, 
which is also a radio emitter, shows associated X- 
ray point source emission from the XMM/Newton- 
COSMOS survey (Brusa et al. 2006). It has 0.5- 
2 keV X-ray emission of 3.74xl0~ 15 ± 0.34xl0" 16 
erg cm s , corresponding to an X-ray lumi- 
nosity of 1.21x10 erg~ 1 s^ 1 at the spectroscopic 
redshift of the lensing galaxy (z — 0.347, Trump 
et al. 2007) . This X-ray emission as well as the op- 
tical spectrum from Magellan (Brusa et al. 2006) 
are consistent with the lensing galaxy being a ra- 
dio galaxy with a central low- luminosity AGN. In 
Table [3] we list the radio counterpart properties 
and show VLA overlays on the ACS images in 
Figure [16] 

5. Simple mass models for the best sys- 
tems 

We want to estimate the Einstein radius of the 
lensing galaxies and determine if there are multi- 
ple images which are not detected, and/or if there 
is a need to take into account for the environment 
to model the system. We use the Lenstool code 
to make a parametrised mass model of the lens 
(Kncib et al. 1996, Julio et al. 2007). The lens- 
ing galaxy is modeled by a Singular Isothermal 
Ellipsoid (SIE). An external shear contribution is 
added when the SIE alone cannot reproduce the 
image configuration. The parameters optimized 
by the code are the coordinates, the orientation, 
the ellipticity and the velocity dispersion of the 
SIE, and the direction and amplitude of the exter- 
nal shear. The (unknown) redshift of the source is 
calculated assuming that the source is at twice the 
comoving distance of the lens. This assumption, 
while clearly incorrect in general, produces a lower 
limit for the estimated velocity dispersion of the 
lensing potential, since the efficiency of lensing is 
greatest for this configuration. 



The Einstein radius is derived using the follow- 
ing equation: 

where D s is the angular distance between the ob- 
server and the source and D[ s is the angular dis- 
tance between the lensing galaxy and the source. 
The constraints used to optimize the mass model 
are the coordinates of the images. The multiplic- 
ity of the systems used in the model is given in 
Table [4] We first model the lens potential with 
a single SIE, and let one or more of the SIE pa- 
rameters vary freely The model uncertainties are 
typically ±0.1 for the ellipticity, and ±10° for the 
orientation, reflecting a possible misalignment be- 
tween the mass and the light (Kochanck 2002). If 
necessary, and if there are enough constraints, an 
external shear contribution can be added. For the 
ring-like objects, we have optimized the velocity 
dispersion of the SIE by hand in order to retrieve 
a perfect Einstein ring. For the purposes of this 
calculation, we assume that the source, the lens 
and the observer are perfectly aligned. Thus we 
have not given any \ 2 value for these cases in Ta- 
ble U The system COSMOS 0013+2249 is not 
included in the table, as it possesses a single vis- 
ible elongated and curved arc, providing too few 
constraints to fit a mass model. 

In Table [4] we display the velocity dispersion 
and the external shear parameters corresponding 
to the best fits. Additional details are given in 
comments in the caption when the result of the 
fitting procedure is not straightforward. 

6. Analysis of the sample and comparison 
to other strong lensing surveys 

In this section we compare the physical parame- 
ters of the foreground in the COSMOS strong lens 
sample with independent galaxy-galaxy strong 
lensing catalogs from the SLACS survey (Treu 
et al. 2005, Bolton et al. 2006) and the CAS- 
TLES database^ (Muhoz et al. 1998). Although 
the SLACS and CASTLES samples were selected 
very differently from the COSMOS strong lens- 
ing sample, they also have HST imaging. Thus a 
rough comparison in terms of the galaxy lensing 



4 http: / /www. cfa.harvard.edu/castles/ 
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Tabic 3: Radio information for the strong lenses. 



Name 


Flux 


S/N 


Flag 


Comment 






(mJy) 










0029+4018 


0.291 


16.710 


1 


likely lensing ; 


galaxy counterpart 


0208+1422 


0.088 


6.960 





likely lensing ; 


galaxy counterpart 


5748+5524 


0.163 


5.680 





likely lensing ; 


galaxy counterpart 


0013+2249 


0.141 


12.450 


1 


likely lensing ; 


galaxy counterpart 


0038+4133 


0.108 


4.740 





elongated, likely lensing galaxy counterpart 


5914+1219 


1.240 


88.520 


1 


strong source, must be a quasar in the galaxy source 










or in the lens 





Note. — Column 1: Lens candidate name. Column 2: Total Flux in mJy. Column 3: Signal to noise ratio of the detection. 
Column 4: Resolved flag (l=resolved radio source, 0=point source). Column 5: Comment. 




Fig. 16. — Radio contours on top of the ACS images. The name on top of each image refers to the radio 
source name in the VLA-COSMOS survey (Schinnerer et al. 2007), and are based on the radio source exact 
coordinates. The contours are defined as follows: 3, 5, 7, 9, 11, 13, 15 and 17cr for COSMOS 0013+2249 
and 0029+4018; 3, 4, 5 and 6cr for COSMOS 0038+4133, 0208+1422 and 5748+5524; 3, 6, 12, 24 and 48ct 
for COSMOS 5914+1219. 
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Tabic 4: Morphological parameters and Einstein radii of the lensing galaxies for the best systems. 



Name 


z/ 


z s 


x 2 


a v 




7,0 7 


More images? 


Multiplicity 










(km s- 1 ) 


(") 




Comments 




0012+2015 


0.41 


0.95 


0.01 


215.3 


0.67 


no 


no 


Double 


0018+3845 


0.71 


1.93 


1.9 


303.1 


1.32 


(0.28,47.5) 


no 


Triple 


0038+4133 


0.89 


2.70 


0.08 


225.3 


0.73 


(0.06,179.1) 


no 


Ring 


0047+5023 


0.85 


2.51 


1.4 


313.0 


1.41 


(0.23,174.1) 


no 


Triple 


0049+5128 


0.33 


0.74 




380.0 


2.09 


no 


no (1,2) 


Ring 


0050+4901 


1.01 


3.34 


3.4 


342.3 


1.69 


(0.23,76.1) 


no (2,3) 


Quad 


0056+1226 


0.44 


1.03 


1.3 


337.4 


1.64 


no 


no 


Double 


0124+5121 


0.84 


2.47 




245.0 


0.86 


no 


no (1) 


Ring 


0211+1139 


0.90 


2.76 


0.01 


466.3 


3.14 


no 


no (2) 


Double 


0216+2955 


0.67 


1.77 


0.06 


348.5 


1.75 


no 


no (2) 


Double 


0227+0451 


0.89 


2.70 


20.0 


428.3 


2.64 


no 


yes (2,3) 


Double 


5857+5949 


0.39 


0.89 


16.2 


398.1 


2.28 


no 


yes (2,3) 


Double 


5914+1219 


1.05 


3.57 


1.2 


338.6 


1.65 


(0.28,68.8) 


yes (2) 


Triple 


5921+0638 


0.45 


1.06 


0.48 


221.0 


0.70 


(0.09,100.5) 


no 


Quad 


5941+3628 


0.90 


2.76 




285.0 


1.17 


no 


no (1) 


Ring 


5947+4752 


0.28 


0.61 




370.0 


1.97 


no 


no (1) 


Ring 



Note. — Column 1: name. Columns 2 and 3: lens and source redshift used for the modeling. Column 4: \ 2 of the fit between 
the best mass model and the data. Columns 5 and 6: velocity dispersion and Einstein radius of the lens. Column 7: External 
shear contribution? if yes, value and orientation of the shear. Column 8: Do more images appear during the modeling? yes or 
no. The labels refer to the comments given hereafter. Column 9: Multiplicity of the images used for the mass modeling. 

Note. — (1) The mass model is created using a source aligned with the lens center, and we scale a v in order to match the 
ring size. (2) The angular distance between the images and the velocity dispersion of the lens makes us think that the lensing 
galaxy is associated with a galaxy group or cluster. (3) The \ 2 value is high, suggesting that this system is not a genuine lens, 
or that a simple mass model is not a good representation of the total lens potential. 
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properties is warranted, at least for the sub-sample 
of best systems. 

6.1. Redshift and effective radius 

Photometric redshifts have been measured by 
Mobasher et al. (2007) for most of the COSMOS 
strong lensing galaxy candidates (see Table [2] and 
0. We first compare the redshift distribution of 
the lensing galaxy candidates found in the COS- 
MOS field (17 best candidates with a photomet- 
ric redshift) to the lensing galaxies in the CAS- 
TLES database (90 strong galaxy-galaxy lenses, 
lens redshift known for 59 of them) and in the 
SLACS survey (19 strong galaxy-galaxy lenses) 
(Fig. 15 1. The (photometric) redshift distribu- 



tion of the COSMOS sample of strong lenses ex- 
tends to significantly higher redshifts than the 
spectroscopic redshift distribution of galaxies in 
the SLACS sample, which is intrinsically limited 
due to the SDSS depth, but is similar to the spec- 
troscopic redshift distribution of galaxies in the 
CASTLES sample. 

The effective radii and arc radius of the lensing 
galaxies (in kpc) are simply obtained by convert- 
ing the effective radii in angular units to physical 
distances, via the photometric redshift of the lens- 
ing galaxy. The SLACS sample includes lenses 
with larger effective lensing radii than the COS- 
MOS sample (see Fig. 17 1. This is probably 



an observational bias, as (for practical reasons) 
we inspect visually the lens candidates in a fixed 
10"xl0", box, and therefore miss any galaxies 
that would generate arcs at larger radii. 

6.2. The absolute magnitude 

The absolute magnitudes of the SLACS lensing 
galaxies are available in Treu et al. (2006). The 
mean value of their V-band magnitude is (My) = 
— 22.71±2;4g. The absolute V-band magnitudes of 
the lensing galaxy candidates of our "best sys- 
tems" sample have a mean value of (My) = 
-21.92±^|. Thus the COSMOS sample reaches 
slightly fainter magnitudes than the SLACS sam- 
ple, despite the fact that it extends to higher lens 
redshifts. 

6.3. The velocity dispersion 

The velocity dispersions of the lensing galaxies 
of the COSMOS sample are derived from simple 



mass models presented in Sect. [5] They are dis- 
played in Table [3] We have compared the veloc- 
ity dispersions of the lensing galaxy candidates of 
the best systems to the velocity dispersions of the 
lensing galaxies in the SLACS and in the CAS- 
TLES samples (also measured using isothermal 
mass potentials, Fig. 18 1. The lensing galaxies in 



the COSMOS sample have a slightly higher veloc- 
ity dispersion than the lensing galaxies in SLACS 
and in CASTLES. The difference with the SLACS 
sample can be due to the fact that we are probing 
in mean larger arc radius (therefore larger Einstein 
radius). In addition, we have estimated the source 
redshift from the lensing galaxy photometric red- 
shifts, that have an uncertainty. This can lead in 
some cases to wrong estimations of the velocity 
dispersion of the main deflector. In particular, we 
notice that two systems have a high velocity dis- 
persion (> 400 km s- 1 ) : COSMOS 0211+1139 
and COSMOS 5914+1219. If real, these values 
suggest that these galaxies are in the center of a 
group, and that we are measuring here both the 
galaxy mass and the group mass. 

6.4. The Einstein radius 

The Einstein radius is directly derived from the 
velocity dispersion of the lensing galaxy, it is there- 
fore another mass estimator for the lens potential. 
It has the advantage of being directly comparable 
to the observable arc radius. The value of the Ein- 
stein radii of the lensing galaxies Re are given in 
Table [4] For the best systems, the Einstein radius 
is expected to be similar to the arc radius, while 
an equality between these two values is expected 
in case of perfect alignment between the observer, 
the source and the lens (Einstein ring). The plot 
in Fig. 19 shows that the two radii are indeed sim- 
ilar for most of the observed systems. Neverthe- 
less we observe some exceptions: some systems 
have a large Einstein radius in comparison to their 
arc radius: COSMOS 0018+3845 (r Qrc =0.40", 
R E =IA4"), 0047+5023 (r Qrc =0.70", R E =1.29") 
and 5914+1219 (r orc =1.86", R E =2.82"). For the 
three systems, there is a strong external shear 
added to the SIE. This may indicates the pres- 
ence of a massive structure in the direction to the 
source or in its vicinity that could have stretched 
the angular separation between the images and/or 
biased strongly the estimation of the velocity dis- 
persion, and hence of the estimation of the Ein- 
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stein radius of the lensing galaxy. 

7. Discussion and Conclusion 

In this paper we present 67 new strong galaxy- 
galaxy lensing candidates found by visual inspec- 
tion of the COSMOS field. From a catalog of 
278,526 galaxies with photometric redshifts (Ca- 
pak et al. 2007), we have selected a sub-catalog 
of 9452 bright galaxies (My <-20 mag), at inter- 
mediate redshift (0.2 < z v h ot < 1.0) and spec- 
trally early type. The selection criteria for this 
sub-catalog were mainly motivated by the pur- 
pose of reducing the sample of galaxies to be vi- 
sually inspected. From the visual inspection of 
ACS/lF8i4i«-band stamp images of 10"xl0", we 
have built a sample of 337 possible lenses with 
sizes up to r arc ~5". This number is reduced to 
67 strong galaxy-galaxy lensing candidates after 
a final cut based on an inspection of the pseudo- 
color images, and the fit and subtraction of the 
lensing galaxy luminosity profile. There are 47 
candidates with a single arclet and 20 candidates 
with multiple images of a background source or 
with a large curved arc, called the "best systems" 
in this paper. We have produced a simple mass 
model (SIE+shear) for each of the best candidates. 
In a few cases, we have also been able to predict 
the existence and direction of an external shear 
contribution. Comparing the velocity dispersions 
to those of lenses in the the SLACS and CAS- 
TLES samples, it is apparent that we are probing a 
slightly higher range of velocity dispersion. Never- 
theless, massive objects in the immediate vicinity 
of the COSMOS lenses, and the rough assumption 
concerning the source redshift, can explain this re- 
sult. We have inferred the Einstein radius for the 
lensing galaxies considering that they are isother- 
mal spheres, and compared these values with the 
arc radius. The derived Einstein radii agree well 
with the observed arc radii, making us confident 
that the sample consists mainly of genuine lens 
systems. 

The "visual inspection" approach is long and 
tedious, and does not ensure a full understand- 
ing of the completeness of the strong lenses sam- 
ple. Nevertheless, our meticulously-obtained large 
sample of strong lenses, even if incomplete, pro- 
vides a valuable catalog to help develop, test 
and improve automatic strong lens-finding algo- 



rithms for similar data-sets. Such work is cur- 
rently in progress (Cabanac et al. 2007, Marshall 
et al. 2007, Seidel & Bartlemann 2007). Also, in 
future work (Faure et al. 2007) the direct study of 
the spatial correlation between large-scale struc- 
ture in the COSMOS field (Scoville et al. 2007, 
Massey et al. 2007) and the strong lenses discov- 
ered here will allow us to measure the impact of 
the total, projected mass distribution along the 
line of sight to a given strong lens, as recently in- 
vestigated for three strong lenses in the AEGIS 
survey (Moustakas et al. 2006). A spectroscopic 
follow-up of the lens candidates is in progress 
at the VLT using the FORSl-Multi-Object Spec- 
troscopy capability. This will allow us to mea- 
sure precisely the redshift of both the lens and the 
source, and to explore the close neighborhood of 
the lens, thus allowing a definite measurement of 
the mass of each strong lensing galaxy. 

Finally, this study gives a lower limit on the 
number of strong lenses expected in future deep 
space surveys such as the proposed JDEM/SNAP 
multicolor weak- lensing survey. Indeed, we expect 
to find at least 10 strong lenses per square degrees. 
Thus for a 10 000 square degree survey, we should 
discover about 100 000 strong lenses thus allowing 
a high precision statistical analysis of the mass 
properties of galaxies and its evolution over time. 
The sample of strong lenses is available at this 
address: |http://cosmosstroir glensing. uni-hd.de/ 
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Fig. 17. — The effective radii of strong lenses 
in the SLACS survey (hashed histogram) and 
the best systems from the COSMOS survey (his- 
togram filled in grey). 



This 2-column preprint was prepared with the AAS IATgX 
macros v5.2. 
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Fig. 18. — The velocity dispersions of strong 
lenses in the SLACS survey (hashed histogram) 
and in the CASTLES database (solid line) com- 
pared to the velocity dispersions of the best sys- 
tems in the COSMOS sample (histogram filled in 
grey). 
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Fig. 19. — The arc radius versus the Einstein ra- 
dius of the best systems. The black stars show 
the ring like candidates. The solid line traces 
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